Color confinement and dual superconductivity in unquenched QCD 
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We report on evidence from lattice simulations that confinement is produced by dual su- 
perconductivity of the vacuum in full QCD as in quenched QCD. Preliminary information 
is obtained on the order of the deconfining phase transition. 

1. Introduction 

Confinement of color and the deconfining phase transition have been observed in numer- 
ical simulations on the lattice. In pure gauge theory, (SU(2) and SU(3)), the transition is 
well understood. (L), the expectation value of the Polyakov line, is a good order param- 
eter and (L), the expectation value of the 't Hooft line, is a good disorder parameter [ |TJ; 
the corresponding symmetries are Zpj, Z^. 

The situation is less clear in the presence of dynamical quarks (full QCD). There the 
symmetry Z^ (Z^) is explicitely broken by the coupling of the quarks. In the chiral limit 
m q = the chiral parameter (ipip) would be a good order parameter for the chiral phase 
transition; however at m q ^ also chiral symmetry is explicitely broken. Moreover it is 
not clear a priori what relation exists between chiral symmetry and confinement. 

For a model QCD with two quarks (Nf = 2) of equal mass m u = = m q the phase 
diagram is schematically represented in Fig. 1. The line of phase transition in the plane 
m q , T is defined by the maxima of susceptibilities [ 0, [3|, among which 

X l = J d 3 x(L(£,0)Lt(0,0)> (1) 
and the susceptibility Xch of the chiral order parameter 

Xch = / d 3 x(^(x, 0)^(0,0)) . (2) 

All of them have a maximum at the same value of T, for a given m q , which defines the line 
in Fig. 1. For m q > 3 GeV, the maxima of xl diverge proportionally to the volume V, 
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indicating a first order transition. At m q 
indications that the transition is second order 



there are theoretical reasons and numerical 
I, El Bl. At intermediate values of m„ 



the susceptibilities studied in [ pj, [| stay finite as V — > oo, and this is interpreted as a 
crossover. 

In the usual approach to phase transitions the free energy (effective lagrangean in the 
language of field theory) depends on the order parameters in a form which is dictated by 
symmetry and by scale analysis. The susceptibilities which keep memory of the order of 
the transition are those of the order parameters. At intermediate quark masses neither 
(L) nor (ipip) are good order parameters. 




Figure 1. Phase diagram for two degenerate flavours, m u = = m q 
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Figure 2. p as a function of N s for T < T c (left-hand side), p as a function of N s for 
T > T c : as iV s — > oo, p — > — oo if the magnetic charge is non-zero and stays constant 
otherwise (right-hand side). 

A good order parameter exists, which is the vacuum expectation value of a magnetically 
charged operator p. (p) ^ signals dual superconductivity of the vacuum [ |^, |6], [7|, |] . 
According to the ideas of ref.'s [ [| p7o[| , dual superconductivity channels the chromoelectric 



3 



-5000 



-15000 



-20000 





I 




= 16 


10000 




•N a 


= 20 






vN s 


= 32 




5 







I 3 



5.0 



5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 6.0 
P 



Figure 3. p peak at different values of N s . The values of m q are chosen such that m q N] 
is kept constant (m q = 0.075 at N s = 16, m q = 0.043 at N s = 20 and m q = 0.01335 at 
N s = 32). 

field acting between a qq pair into an Abrikosov flux tube, with energy proportional to 
the length, V = or. 

In quenched QCD [ || [7], || a detailed analysis of (p) shows that the vacuum is indeed 
a dual superconductor for T < T c ((p) ^ 0), and makes a transition to normal at T c . 
Above T c , (p) = 0. As T -> T c __, (/i) ~ r 5 0(a/£, CM), with r = (1 - T/T c ), £ ~ r"" 
the physical correlation length and N s the linear spatial extension of the lattice. At the 
critical point a/^ C 1 and can be put to zero, while the variable £/iV s can be traded with 
N]/ v t. As a consequence the scaling law (/i) = r 5 (f)(N^ u T) holds. A better variable than 
(/i) is p = log(/i), which is a susceptibility. In terms of p, 

(//) = exp U '%{P')dP^ (3) 

The scaling law p/N}/ u = f(rN}/ u ) allows one to control the infinite volume limit N s — > 
oo, which is necessary to have a phase transition, and to extract u, T c , 5. For quenched 
QCD, (p) is equally good as (L) and is numerically coincident with it [ 0]. 

A similar analysis can be repeated for full QCD [ [Q]|. (p) is defined as in quenched and 
defines the same symmetry. Simulations have been done on N t x lattices with N t = 4 
and N s = 12, 16, 20, 32. The time extension Nt determines the transition temperature, 
the space extension N s is used to perform the finite size scaling analysis of the infinite 
volume limit. The first result is that for T < T c , (p) ^ as N s — > oo. Indeed in that 
range of temperature, p is N s independent (see Fig. 2). For T > T c instead p ~ —kN s + c 
{k > 0) as N s — > oo (see Fig. 2), i.e., by Eq. (^), (p) is strictly zero in the infinite volume 
limit. 

Around T c , again, 

(p) = T S <p(- c ^,m q Nl) . (4) 
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A new scale, m q , appears which was absent in the quenched case. However the exponent 
7 is known (7 ~ 2.49 [ || |3|) and simulations can be made choosing m q and N s such that 
m q N] in Eq. is kept fixed. Then, neglecting a/£ again, the scaling law becomes 

<A*> = r 5 4> (0, ^, const.) , p/N^ = f(rN^) . (5) 

The behaviour of p around T c is shown in Fig. 3 for different values of N s : it has a peak 
at T c , Eq. (|5|) implies in particular that the values of p at the peak (r = 0) p, scale 
as p oc N^ u , whence v can be determined. The result is shown in Fig. 4: contrary to 
the common belief, the transition is first order. This result is still preliminary and will 
be cross-checked by more simulations, and by the use of improved actions. A scenario 
in which confinement is controled by one single order parameter (p) and by one and the 
same symmetry pattern (dual superconductivity) is appealing and meets the requirements 
of N c — ► 00 arguments. 



REFERENCES 

[1] L. Del Debbio, A. Di Giacomo, B. Lucini, Phys. Lett. B500, 326 (2001). 

[2] F. Karsch and E. Laermann, Phys. Rev. D 50, 6954 (1994). 

[3] S. Aoki et al. (JLQCD collaboration), Phys. Rev. D 57, 3910 (1998). 

[4] R. Pisarski and R. Wilczek, Phys. Rev. D 29, 338 (1984). 

[5] L. Del Debbio, A. Di Giacomo, G. Paffuti, Phys. Lett. B 349, 513 (1995). 

[6] A. Di Giacomo, B. Lucini, L. Montesi, G. Paffuti, Phys. Rev. D 61, 034503 (2000). 

[7] A. Di Giacomo, B. Lucini, L. Montesi, G. Paffuti, Phys. Rev. D 61, 034504 (2000). 

[8] J.M. Carmona, M. D'Elia, A. Di Giacomo, B. Lucini, G. Paffuti, D 64, 114507 (2002). 

[9] G. 't Hooft, in High Energy Physics, ed. A. Zichichi (EPS International Conference, 

Palermo 1975). 
[10] S. Mandelstam, Phys. Rep. 23C, 245 (1976). 



5 



[11] J.M. Carmona, M. D'Elia, L. Del Debbio, A. Di Giacomo, B. Lucini, G. Paffuti, Phys. 
Rev. D 66, 011503 (2002). 



